Quantum metrology: why entanglement? 
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We show why and when entanglement is needed for quantum-enhanced precision measurements, 
and which type of entanglement is useful. We give a simple, intuitive construction that shows how 
entanglement transforms parallel estimation strategies into sequential ones of same precision. We 
employ this argument to generalize conventional quantum metrology, to identify a class of noise 
whose effects can be easily managed, and to treat the case of indistinguishable probes (such as 
interferometry with light). 



To measure a parameter (p that identifies some phys- 
ical transformation t/^, we can sample the transforma- 
tion with a probe system and then measure the change 
in the state of the probe. To reduce the statistical er- 
ror, we repeat the sampling N times and average the 
results. Quantum metrology encompasses a large class 
of estimation procedures where quantum effects such as 
entanglement and squeezing are used to enhance the pre- 
cision of such measurement over what would be possible 
with classical resources JU |j. Such techniques are by 
now well established |q-pO| and have produced impres- 
sive results pl|-^. Yet, the reason why entanglement 
is necessary has not been settled. Here we present a 
simple, intuitive argument on why and when entangle- 
ment is necessary and, additionally, use it to derive some 
new results in quantum metrology. Our argument shows 
how entanglement can convert a parallel-sampling strat- 
egy (where N probes sample the system in parallel) to a 
sequential strategy, which is known to achieve high pre- 
cision, and viceversa. While the equivalence of parallel 
and sequential strategy was already shown in Ref. |T^ 
using the theory of tensor networks, our construction is 
much simpler, as it is based on a trivial algebraic identity: 
this simplicity is key to give an intuition on what quan- 
tum metrology means, to obtain some new results, and 
to reobtain more easily some known ones: e.g., we show 
how one can easily generalize the conventional quantum 
metrology framework; which kinds of errors are simple to 
manage; how the theory accommodates indistinguishable 
probes (e.g. bosons for light interferometry); which types 
of entanglement are useful; how non- asymptotic results 
arise in quantum metrology; etc. 

In the framework of quantum metrology we are given 
a "black box" that performs a unitary operation U^p ~ 
e"^^, where (p is the parameter to be estimated and H 
is its generator of translations. What is the best preci- 
sion with which we can determine ip7 Call |0) and |1) 
the states that correspond to the minimum and maxi- 
mum eigenvalues of H. The existence of such states is 
guaranteed only in finite-dimensional Hilbert spaces, and 
we will restrict to this case here (see Refs. p^, p3] for 
some infinite-dimensional results) . If we sample the sys- 
tem sequentially N times with a single probe (Fig. IlK), 



we can perform a measurement of ip as follows: (a) pre- 
pare the probe in the state |-|-) or |— ), |±) = |0) ± |1) 
(we will neglect state normalizations); (b) evolve it us- 
ing a sequence of A'' unitaries as U^\±) = |0) ± e*^'''|l); 
(c) project the final state onto the initial one: the prob- 
ability that they coincide is p{ip) = cos{Nip/2) (Ram- 
sey interferometry corresponds to the case A^ = 1 Q). 
By repeating the above procedure v times, the experi- 
mental probability p('/3), whence (p, are inferred with an 
error that scales as l/{Ny^) M M where the factor 
Y^ comes from the central limit theorem and the factor 
N comes from the fact that the sequential procedure is 
equivalent to measuring the parameter Nip: an error A 
on Nip corresponds to an error A/N on tp. In contrast, 
if we prepare the Nv probes in the state |+) or |— ), use 
them to sample the system in parallel and average the re- 
sults (Fig. Illb), from the central limit theorem if follows 
that we achieve an error that asymptotically scales as 
I/VNv- the "standard quantum limit" (SQL). Surpris- 
ingly, if the N probes are prepared in an entangled state 
(Fig. yjc), the quantum Cramer- Rao bound can be used 
to show that one can asymptotically achieve the same 
1/A^ limit [Q-^ of the sequential strategy- the "Heisen- 
berg bound". A typical quantum metrology experiment: 
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(a) prepare N probes in the entangled state |0) 

(b) evolve it by applying U^ to each probe in parallel: 
U»^{\0)^ + |1)^) = |0)^ + e^^^ll)^; (c) project on the 
initial state as above to estimate the parameter Np: as 
before, the probability is p{p) = cos{Np/2), whence p 
can be estimated with an error that scales as l/(iV-y/i/) 
as in the sequential case above. While this is well known, 
the reason why entanglement helps was not explored in 
depth: typically one would say that this is the effect of 
the "global" nature of the entangled state that cannot 
be factored into separate contributions for each probe. 
Here we give a simple construction that shows how one 
can convert a parallel-entangled strategy into a sequen- 
tial one and viceversa, proving that this is a reason why 
they achieve the same error scaling. We show that clas- 
sical correlation among the probes is insufficient. 

The argument: — We start detailing our construction 
for A^ = 2 and then extend it to arbitrary A^ by induction. 
Given an operator C = ^- ■ Cij\i){j\ {{\i)} a basis), we 
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FIG. 1: Possible strategies for estimating the parameter ip 
that a unitary transformation U^p (black boxes) encodes into 
a probe system (triangles), (a) Sequential strategies: a sin- 
gle probe samples the A*' boxes sequentially, effectively mea- 
suring Nip and achieving an error oc A''"^ on the estimation 
of if. (b) Classical (unentangled) parallel strategies: we use 
TV probes and average the measurement results: the central 
limit theorem tells us that the error on ip scales as N~^'^. 
(c) Entangled strategy (quantum metrology): the A'^ probes 
are prepared in a joint entangled state, but sample the black 
boxes separately: the error scales as N^^ . 



can define a state \C) = ^- ■ Cij\i)\j). It can be used ||26 
for the following identity 



A^B\C) = J:^m^k^Bl,C.,\k}\^) = \ACB^) 



(1) 



where ^ indicates the transpose in the \i) basis. Consider 
the quantum metrology protocol (Fig. ft) for N ~ 2, and 
apply this identity twice: 



{u^®u^,)\t) ^W^u"^,) ^ {u^u'!^, 
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where |1) = |00) + |11), and we consider the general case 
where ip and if' may differ. Choose \i) as the eigenbasis of 
H, so U^ = e"^-^ is diagonal, and C/J ~ U^. Moreover, 
we can exploit the entanglement of |1) to measure the 
second system in the {|+), |— )} basis (we consider only 
the subspace spanned by |0) and |1)) and project the first 

system in the same state, since |00)-|-|II) = | ++) + | ). 

This means that after the measurement, the evolution 
in Eq. (0) can be written as C/;^L/(^'|±), where the initial 
state is |-|-) or |— ) with probability 1/2, depending on the 
measurement outcome. This argument shows how the 
parallel entangled strategy {U^p (S) Utpi)\t) can be reduced 
to the sequential one U^U^'\±) for N = 2, see Fig. 0. 
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Then, using the eigenbasis of H and projecting 
the second probe into the |±) basis: 
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FIG. 2: Argument to transform a parallel entangled strategy 
into a sequential one for N — 2, see Eq. (B). 

Why is entanglement necessary? We must require that 
perfect correlation is present in two complementary ba- 



sis, the |0, 1) basis (to ensure that 11,^ is diagonal) and 
the |±) basis (to ensure that the measurement on one 
probe projects the other in an appropriate initial state). 
An intuitive definition of entanglement is "correlation on 
a property that does not (cannot) yet exist" , namely, 
a perfect correlation on complementary properties such 
as the two bases above. In fact, suppose we request 
that there is correlation only on the |0, 1) basis using 
the mixed state |00)(00| + |11)(11|), that can be written 
as an equally weighted mixture of |1) = |00) + |11) and 
\<Jz) = |00) — 1 11), each with probability 1/2. We can use 
the construction of Eq. (pi) obtaining a final state of the 
two probes given by [U^U^ ® 1)|1) or [11^(7^11^ (g) 1)|1) 
each with probability 1/2. After the projection of the 
second probe into the |±) basis, these become an equally 
weighted mixture of |0)±e*2'^|l) and |0)Te*^'^|l), namely 
the state |0)(0| -I- |1)(1|, independent of ip. In contrast, 
suppose we request correlation only on the |±) basis us- 
ing the state | ++)(++ | + | )( |, namely an 

equally weighted mixture of |1) = | + +) + | ) and 

\<Jx) = 1+ +> - I- -> (where a, = |0)(1| + |1)(0| in the 
|0,1) basis). Again using Eq. (g), this is equivalent to 
having a final state of the two probes as {U^U^ (E) 1)|1) 
or {U^<JxUtp (g) 1)|1) each with probability 1/2, which, 
after the projection, become an equally weighted mix- 
ture of the four states |0) ± e^^'^jl) and |0) ± |1), again 
|0)(0| + |1)(1| independent of ip. [Here the SQL is re- 
covered if one does not average over the projective- 
measurement result: that procedure is equivalent to the 
classical parallel strategy of Fig. |l|b.] 

We have shown that for A^ = 2 we can easily convert 
a parallel entangled strategy into a sequential one, prov- 
ing that the two have the same estimation precision. We 
now extend this to arbitrary A'^ by induction: suppose 
(induction hypothesis) that the parallel entangled strat- 
egy 07=1 ^vjil^)^ + |1)^) has the same estimation pre- 
cision as the sequential strategy Jlfci ^vj\^)j we need 
to prove that the same applies for iV -t- 1. We can write 
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where |0) = |00), |1) = |11), and 
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since C/;^ 's are diagonal in the |0, 1) basis. It follows 
that the subspace spanned by |0) = |00) and |I) = |11) 
is invariant for application of C/g, which, when restricted 



to such subspace, can be expressed as the 2x2 matrix along the hnes of Eq. (0) as 



VO e^^f"~^'^"+^^) ~ ^vnU<pk+i ■ (5) 



Restricting ourselves to such invariant subspace, we can 
use the induction hypothesis to conclude that the last 
line of (@) will have the same estimation precision as 
Yij^i^ U^^Ue\±) = rij-V f^vji)' which proves the in- 
ductive step, see Fig. y. [Note that in this proof it is 
not necessary to require that the matrix Ug is diagonal, 
but only that it is invariant for the |0), |1) subspace: an 
anti-diagonal tensor product would work.] 
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FIG. 3: Extension of the argument of Fig. H to arbitrary N. 



Consequences: — We now explore some of the conse- 
quences that can be derived from our argument. 

(a) Useful entanglement: Up to now we have used 
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as entangled initial state the GHZ, (|0) 
optimizes the quantum Cramer-Rao bound. Our con- 
struction shows that entanglement of the type I^at) = 
|0)^ -f e''^|l)^ (with real A) is the only useful one (with 
the appropriate amendments when probes are indistin- 
guishable, see below). In fact, suppose we use an arbi- 
trary entangled state \E) in Eq. (0), it becomes (e"^'^ (g) 
e''<'")\E) = {e"f"Ee"^" ® l)\t), where we want the 
right hand term to correspond to a sequential strategy 
^2ivHQQ^ ± e*^|l)) (the sequential strategy works iff the 
initial state is of the form 1 0) ± e*'^ 1 1) even though for sim- 
plicity we considered only the cases |±) above). Namely, 
we require that e'^^" Ee'vH\±) ^ e^'f"{\0) ± e'^|l)), 
which implies that E == |0)(0| + e''^|l)(l|, since H is 
diagonal in the |0,1) basis. Namely, we need a state 
\E) — 1^2)- This argument, and the fact that achieving 
the quantum Cramer-Rao bound requires a superposi- 
tion of global maximum and minimum eigenstates [pi, 
shows that quantum metrology requires entanglement of 
the form \^n)- 

(b) Noise resilience: using our framework, we can de- 
tail a class of noise that is simple to handle (in general, 
the treatment of quantum metrology in the presence of 
noise is highly nontrivial p5|-p^, p7|-p9t). Consider the 
entangled parallel strategy for iV = 2 and suppose that 
the two probes are subject to arbitrary noise maps de- 
scribed by Kraus operators Ak and B^ respectively. We 
can still convert a parallel strategy into a sequential one 



Y,{Ak® B,)\l){t\{Al® b]) ^Y.\'^>^Bj){AlB]\ 

jk jk 

= 5](A,Bj®i)|i)(i|(i?;4®i), (6) 

jk 

i.e. a noise map only on the first probe if ^ BjB^, ~ 1 
(a unital map, as when Bj ex (7^, the Pauli matrices), 
since then J2jkB*AlAkBj = 1, with B* = (B'^)t. So 
the first step of the induction can be applied also to uni- 
tal maps in addition to the unitaries we discussed above. 
However, the extension to A^ > 2 fails in general, as there 
is typically no invariant subspace to construct a reduced 
map as the one in (|5|): one must require that A^ ^ Bj is 
a diagonal or an anti-diagonal matrix for all k,j. An ex- 
ample of the latter is the bit flip with phase noise where 
the Kraus operators are Aq oc ax and Ai oc ay. An ex- 
ample of the former (in addition to the noiseless case) is 
phase flip noise (dephasing) with Kraus operators Aq (x 1 
and Ai oc az, which was already identified in |ll| as a 
noise model whose parallel and sequential performance 
match. It was known for a long time [n7| that an arbi- 
trary small amount of dephasing would ruin any advan- 
tage for frequency measurement [|l^, [l^-[2q| whereas the 
same identical noise does not affect phase estimation as 
much [^, gO| . Our construction explains this curious per- 
formance: since the entangled parallel strategy is (for this 
noise model) equivalent to the sequential one, it means 
that the entangled strategy's performance is equivalent 
to one that takes TV time as long: namely, its noise sensi- 
tivity is TV times faster |0|. While this is not a problem 
for estimating a phase (which can be done in arbitrar- 
ily short time), it is a problem for estimating frequency 
which requires a measurement duration that voids any 
advantage from entanglement [0, O. EH]. The calcula- 
tion of the quantum Cramer-Rao bound for frequency- 
estimation on a single qubit subject to dephasing (the 
optimal sequential strategy) is given in Ref. [Q , whence 
one can see it matches with the optimal entangled strat- 
egy (although that was not observed there). Our con- 
struction thus identifies some noise models whose effect 
is equivalent for the sequential and parallel strategies: the 
unital maps (in the |0), |1) invariant subspace) where all 
Kraus operators are either diagonal or anti-diagonal, as 
the examples above. It is then simple using our construc- 
tion to find the states that optimize the precision in the 
presence of these classes of noise. 

(c) Indistinguishable probes: the above discussion fo- 
cused on the case of distinguishable probes, where a ten- 
sor product structure (TPS) is clear: each probe is a sys- 
tem with its own Hilbert space. If the probes are indistin- 
guishable particles (e.g. in optical interferometry, where 
typically one considers single-photon probes) , such struc- 
ture is absent, and one typically resorts to the Fock space 
that enforces the symmetry properties required (symme- 



try for exchange of Bosons and anti-symmetry for ex- 
change of Fermions). In the absence of a TPS, the very 
concept of entanglement is ih defined [|o[ |3^, and it is 
known that entanglement is not necessary to beat the 
SQL in interferometry. For example, the unentangled 
"NO" state |iV) 4- 10) , where |iV) is a Fock state of N pho- 
tons and |0) is the vacuum state (we restrict ourselves 
to a finite-dimensional subspace of the radiation Hilbert 
space) can achieve a phase sensitivity beyond the SQL: its 
evolution through e"^° ° (a being the annihilation opera- 
tor of the mode) yields |0) -|-e'^'^|A^), as in the sequential 
or entangled strategies of Figs. ||a,c. Because such a state 
cannot be attained in the presence of super-selection rules 
p2| and it only samples an absolute phase, in quantum- 
enhanced interferometry one typically considers the two 
mode "NOON" state g |iV,0) + |0,iV), which samples 
a relative phase between the two modes a, b of the in- 
terferometer through the evolution e*"^^" ""'' ''■'. While 
clearly |A^, 0) + |0,-/V) has mode entanglement, the state 
|A^) -I- |0) has no entanglement. Yet, they can both be 
easily cast in the framework we have used up to now, 
by recalling that for Boson probes, one needs to con- 
sider only the subspace of the global Hilbert space that 
is symmetric for exchange of the probes, which can be 
done by moving to the Fock space Q: in each mode 
we have that the tensor product of minimum and max- 
imum eigenstates |0)^ and |1)^ is symmetrized as |0) 
and I TV) , respectively. Analogously, in the two mode case, 
the symmetrized minimum and maximum eigenstates for 
the two-mode Hamiltonian a^a — b^'b restricted to the A^- 
photon subspace is |0,A^) and |iV, 0), respectively (the 
symmetrization must be performed in each mode sepa- 
rately, since photons in different modes are distinguish- 
able). It is hence clear that NOON state interferometry 
falls within the framework described here. In particular, 
our construction shows that NOON state interferometry 
is indeed equivalent to multipass interferometry, as was 
pointed out previously |lj, |25| . The same trick of consid- 
ering the (anti-symmetric) Fock space can also be applied 
to the Fermionic case, although the Pauli exclusion prin- 
ciple makes Fermionic probes unattractive. [One may 
wonder whether the conventional theory of distinguish- 
able probes is ever applicable, since clearly all probes will 
be ultimately be composed of indistinguishable particles. 
We remind that indistinguishable particles are subject 
to specific symmetries only regarding their global state. 
Typically one uses only certain degrees of freedom of the 
probe system for measurement purposes and the remain- 
ing degrees of freedom (spin, position, etc.) can be used 
to enforce the required symmetry.] 

(d) Extension of quantum metrology: Conventionally, 
in quantum metrology the transformation that acts on 



dependent of (p. In this case the sequential strategy 
cannot be a simple iteration of U' as that will not 
typically increase the parameter ip hy N times (e.g., 
when W — t, V = ax ^c find that U^ ex 1), but 
^iNvH ^ (M^t[/^i/t)W. This suggests that a general- 
ized entangled parallel strategy consists in the evolution 
(W^1'[/^V^t)®W(|o)^ + |1)^), which corresponds to the se- 
quential strategy {W''U'^V'')^\±). 

Is our framework general? The only hypothesis made 
so far is that there exists eigenvectors |0) and |1) for H 
connected to the maximum and minimum eigenvalues. 
This requirement, certainly true in finite-dimensional 
Hilbert spaces, implies that the states |-|-) and |— ) (for 
the sequential strategy) and the GHZ state |0)^ 4- |1)^ 
(for the parallel strategy) arc ones with maximum spread 
AH, that saturate the quantum Cramer-Rao bound |3[- 
g. Our framework is thus general for finite-dimensional 
systems, and allows one to construct new quantum 
metrology protocols [g|: all that is needed is a GHZ- 
type state with a superposition of maximum and mini- 
mum eigenstates of H to achieve a quantum metrology 
protocol to estimate the parameter (p of which H is the 
generator of translations. 

(e) The non-asymptotic case: all results given up to 
now (and most of the literature) refer to the asymptotic 
case where the number of repetitions i^ — > oo. How- 
ever, thanks to our construction, it is possible to adapt 
the non-asymptotic results |g, |4| ^, 34 obtained for 
the sequential estimation (by adapting the Kitaev phase- 
estimation algorithms) also to parallel entangled strate- 
gies. This equivalence was already pointed out in Ref. [g: 
all these strategies employ multiple rounds of sequen- 
tial iterations or, equivalently, multiple rounds of NOON 
states with different N. A completely different strategy 



is given in |22 



the probes takes the form U^ 



ai-vH ^ Our construc- 



tion makes it easy to generalize this to a more gen- 
eral U' = We^'^^V, where W and V are unitarics in- 



(f) Bandwidth-noise tradeoff: our construction clarifies 
that the true power of entanglement in quantum metrol- 
ogy is in allowing one to achieve the same vA precision 
enhancement of the sequential strategy, while keeping the 
same A"-fold sampling-time gain of the parallel strategy 
over the sequential one |35|] : high precision in short time. 

Conclusions: — We presented a simple construction 
that shows why entanglement is necessary in quantum 
metrology: to transform a parallel strategy into a se- 
quential one, we need perfect correlation in complemen- 
tary bases: both in the eigenbasis |0, 1) of the generator 
H and in its complementary |±) basis. Such correlation 
requires entanglement. We also presented some applica- 
tions: our construction easily derives some new and some 
known results in quantum metrology. 

Acknowledgments: — I acknowledge useful discussions 
with V. Giovannetti and support by the Royal Society for 
attending the International Scientific Seminar on Sources 
and signatures of quantum enhancements in technology. 



[1] V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 

1330 (2004). 
[2] V. Giovannetti, S. Lloyd, and L. Maccone, Nature Phot. 

5, 222 (2011). 
[3] S. L. Braunstein, C. M. Caves and G. J. Milburn, Ann. 

Phys. 247, 135 (1996). 
[4] S. L. Braunstein, C. M. Caves, Phys. Rev. Lett. 72, 3439 

(1994). 
[5] V. Giovannetti, S. Lloyd, L. Maccone, Quantum metrol- 
ogy, Phys. Rev. Lett. 96, 010401 (2006). 
[6] D.W. Berry,B.L. Higgins, S.D. Bartlett, M.W. MitcheU, 
G.J. Pryde, and H. M. Wiseman, Phys. Rev. A 80, 
052114 (2009). 
[7] S. Boixo, A. Datta, M. J. Davis, A. Shaji, A.B. Tacla, 

CM. Caves, Phys. Rev. A 80, 032103 (2009). 
[8] H. Lee, P. Kok , J.P. Bowling, J. Mod. Opt. 49, 2325 

(2002). 
[9] S. Boixo, A. Datta, S.T. Flammia, A. Shaji, E. Bagan, 
CM. Caves, Phys. Rev. A 77, 012317 (2008). 
[10] D.J. Wineland , J.J. BoUinger, W.M. Itano, F.L. Moore, 

Phys. Rev. A 46, R6797 (1992). 
[11] S. Boixo, C. Heunen, Phys. Rev. Lett. 108, 120402 

(2012). 
[12] C Gross, J. Phys. B: At. Mol. Opt. Phys. 45, 103001 

(2012). 
[13] CM. Geremia, J.K. Stockton, H. Mabuchi, Phys. Rev. 

Lett. 94, 203002 (2005). 
[14] B.L. Higgins, D.W. Berry, S.D. Bartlett, H.M. Wiseman, 

G.J. Pryde, Nature 450, 393 (2007). 
[15] B.M. Escher, R.L. de Matos Filho, L. Davidovich, Nature 

Phys. 7, 406 (2011). 
[16] M. Kacprowicz, R. Demkowicz-Dobrzaiiski, W. 
Wasilewski, K. Banaszek, LA. Walmsley, Nature 
Phot. 4, 357 (2010). 



[17 

[is; 

[19 

[20 
[21 
[22' 
[23| 

[24 

[25' 
[26; 

[27 



[28; 

[29 

[30 
[31 

[32; 

[33 

[34 

[35; 



S.F. Huelga, , et al. Phys. Rev. Lett. 79, 3865 (1997). 
R. Demkowicz-Dobrzaiiski, J. Kolodynski, M. Guta, Na- 
ture Comm. 3, 1063 (2012). 

D. Ulam-Orgikh M. Kitagawa, Phys. Rev. A, 64, 052106 
(2001). 

A. Shaji, CM. Caves, Phys. Rev. A 76, 032111 (2007). 
M.W. MitcheU, Proc. SPIE 5893, 589310 (2005). 
M. Hayashi, Progress in Informatics 8, 81 (2011). 
V. Giovannetti, S. Lloyd, L. Maccone, Phys. Rev. Lett. 
108, 260405 (2012). 

V. Giovannetti, L. Maccone, Phys. Rev. Lett. 108, 210404 
(2012). 

A. Luis, Phys. Rev. A 65, 025802 (2002). 
G.M. D'Ariano, P. Lo Presti, M.F. Sacchi, Phys. Lett. A 
272, 32 (2000). 

U. Dorner, R. Demkowicz-Dobrzahski, B.J. Smith, J.S. 
Lundeen, W. Wasilewski, K. Banaszek, LA. Walmsley, 
Phys. Rev. Lett. 102, 040403 (2009). 
T.W. Lee, S.D. Huver , H. Lee, L. Kaplan, S.B. Mc- 
Cracken, C Min, D.B. Uskov, C.F. Wildfeuer, G. Veronis 
& J.P. Dowling, Phys. Rev. A 80, 063803 (2009). 
S. Knysh, V.N. Smelyanskiy, G.A. Durkin, Phys. Rev. A 
83, 021804(R) (2011). 

P. Zanardi, Phys. Rev. Lett. 87, 077901 (2001). 
P. Zanardi, D.A. Lidar, S. Lloyd, Phys Rev Lett 92, 
060402 (2004). 

S.D. Bartlett, T. Rudolph, R.W. Spekkens, Rev. Mod. 
Phys. 79, 555 (2007). 

P. Teller, An Interpretive introduction to quantum field 
theory, (Princeton Univ. Press, 1995). 
B.L. Higgins, D.W. Berry, S.D. Bartlett, M.W. MitcheU, 
H.M. Wiseman, G.J. Pryde, New J. Phys. 11, 073023 
(2009). 

S. Braunstein, G. Milburn, (separate) private communi- 
cations. 



